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ABSTRACT 

High-velocity features (HVF), especially of Ca II, are frequently seen in Type la supernovae 
observed prior to B -band maximum (Bmax). These HVF start at more than 25,000 kms -1 
in the days after first light, and slow to about 18, 000 km s -1 near Bmax. To recreate the 
Ca II near-infrared triplet (CaNIR) HVF in SN 2011fe, we consider the interaction between a 
Type la supernova and a compact circumstellar shell, employing a hydrodynamic 1—D simulation 
using FLASH. We generate synthetic spectra from the hydrodynamic results using syn++. We 
show that the CaNIR HVF and its velocity evolution is better explained by a supernova model 
interacting with a shell than a model without a shell, and briefly discuss the implications for 
progenitor models. 

Subject headings: supernovae: general — supernovae: individual (SN 2011fe) — line: formation — line: 
profiles 


1. Introduction 

Type la supernovae (SN la) provide a fun¬ 
damental tool for our understanding of the his¬ 
tory of the universe. SN la are ‘standardizable 
candles’ used to explore the expansion of the 
universe as well as the chemical enrichment of 
galaxies dRiess et al.| 1998 : Perlmutter et al. I il999t 


Tsuiimoto fe Shigevamal 120121. The configura¬ 
tion of the progenitor system and the cause of 
the explosion remain elusive. Observations of 
SN 2011fe within the first day after the explo¬ 
sion have shown that any optically thick material 
was within O.I R^ of the exploding star prior to 


was witmn u. utfr, ot tne e xploding star prior to 
the e xplosion ( Nugent et al.l[2oTH |Piro fc Na. kail 
20141) . Spectroscopy of SN la in the first days 
and weeks after the explosion also reveal high 


velocity f eatures (HVF) in Ca II, Si II and 
other io ns (Hatano et alll999t Parrent et al.l20l3 


iMarion et all 2013 ). The HVF has a velocity 

;> 25,000 km s -1 at 15-18 days before B-band 
maximum (Bmax) and slows to a plateau of about 
18, 000 by Bmax. The photosphere also starts at 
high velocity and moves to lower velocities over the 
same interval, but the HVF consistently remain > 
7, 000 kms -1 fa ster than the photospheric velocity 
featu r es (PVF) (IM arion et al.ll2013 ; Maguire et al 


2014 : Silverman et al. 20151) . Ca II HVF ap¬ 


pear i n over 90% of normal SN la ( Mazzali _eL_aL 
2014: Maguire et al.1 2014 


2005: Childress et al. 

2014; 

Silverman et al. 2015 

) and 

larization 

Wang et al. 2003}) 


material has a high covering factor and is asym¬ 
metric. 
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The source of the high-velocity material will 
give insight to the nature of the progenitor system 
or mechanism by which the explosion is initiated; 
understanding both of these is necessary to control 
the systematics in the use of SN la as cosmologi¬ 
cal probes. Previous suggestions for the source of 
the high-velocity ma terial include plum es of par¬ 
tially burned ejecta ( Wang et al • 2003), buoyant 
bubb les due to a gravitatio nally confined detona¬ 
tion ( Kasen fe Plewal[2005lh or interaction with a 
circumstellar medium (CSM) with a total mass 
of high-velocit y material of about 0 . 02 Mp, and so¬ 
lar abundance ( Gerardv et alJ[2004 ; Quimbv et al 


l2006l h These models consider HVF velocity at the 
time near or after Bmax, when the HVF are in 
the asymptotic phase and are fading. It is neces¬ 
sary to explain the full evolution of the velocity 
and strength throughout the period the HVF is 
detectable in order to understand the source. 


In this letter, we consider the interaction be¬ 
tween a circumstellar shell (CSS) of mass 0.005Mq 
located within 0.3 Rq of the center of the explo¬ 
sion and compare this model to a supernova with 
no CSS sources of the HVF. We concentrate on the 
Ca II near-infrared triplet (CaNIR) in SN 201 lfe 
between —16.5d and +2.9d. In (J21 w e describe the 
hydrodynamic methods used, the SN model ( £12.11) . 
and the CSS (' £12.21) . In £j3]we discuss the methods 
used for generating synthetic spectra from the hy¬ 
drodynamic results. We present the results in 
and summarize our conclusions in £j5] 


EOS and the desire to have similar end times for 
each model. 


2.1. Explosion Model 


T he exp losion models are those of Gamezo et 
al. 020051 ) delayed-detonation models b and c. 
We maintain the resolution of the original model, 
which has a maximum radius of 5.4 x 10 8 cm re¬ 
solved into 128 3 zones (4.2 x 10 s cm per zone). We 
use spherical averaging to reduce the model from 
3-D to 1 D. The model provides C, O, Mg group, 
Si group, and Fe group abundances for each zone. 
Exterior to the ejecta we apply a CSM with den¬ 
sity 10 -9 g cm -3 and temperature of 20, 000 K. 

2.2. Circumstellar shell model 

We have explored a range in masses and den¬ 
sity distributions for the shell. These results 
that w ill be reporte d i n a forthcoming paper 
( Mulligan fe Wheeler in preparation!) . For this 
work, we use a shell with inner radius 7 x 10 8 cm 
and outer radius 2 x 10 10 cm (0.3Rq), a total mass 
of 0.005 M© with a density profile that decreases 
linearly outward, and a temperature of 30, 000 K. 
Exterior to the shell, the CSM parameters are 
applied as described in the previous subsection. 
For the hydrodynamic simulation, we assume a 
hydrogenic shell with solar abundances. Total He 
abundance in the shell is taken to be solar, but we 
enhance the 3 He by 10% to act as a tracer. 


2. Hydrodynamics 


3. Spectral synthesis 


We use FLASH 4.1 ( jFrvxell et al . 2000 ) with 
multi pole gravity, 1 D spherical geometry , the 
Hel mholtz equation of state (EOS) (jTimmes & Arnett 
1999h for the supernova-only model and the 


gamma law EOS for the model including the CSS. 
The simulation volume has the supernova at the 
center and a radius of 10 12 cm, with a minimum 
resolution of 4.2 x 10 6 cm per zone. The multi¬ 
species unit is used to track H, 3 He, 4 He, 12 C, 
16 0, 24 Mg, 32 Si, and 56 Ni. The model is evolved 
to about 50s after the explosion at which time 
the the shock has fully propagated through the 
shell (when it is included) and both the ejecta and 
shell are in or near the free expansion phase. The 
choice of end point and EOS for a given model 
is governed by the density floor in the Helmholtz 


We use the result of the hydrodynamic simu¬ 
lation to generate an optical depth profile as a 
function of velocity. The line optical depth for 
an individual zone is computed using the Sobolev 
approximation then binned with zones having a 
similar velocity. We assume excitation tempera¬ 
ture is constant over the line-generating regions 
in the ejecta and shell. We define scaling factors 
C E and C s to account for the fraction of Ca II 
relative to all Si group elements in the ejecta and 
Ca II abundance by mass within the shell, respec¬ 
tively. 


We use syn++ ( Thomas et al . 201 ill , modified 
to allow an arbitrary optical depth as a function of 
velocity, to generate the synthetic spectra. When 
generating the spectra, the time after explosion, 
photospheric temperature and velocity, and the 
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abundance scale factors C E and C s are free pa¬ 
rameters. The photospheric velocity and tempera¬ 
ture set the minimum velocity of the photospheric 
line-generating region and the shape of the con¬ 
tinuum, respectively. We adopt a lower limit to 
the temperature of the photosphere of 8, 000 K. 
The effects of the excitation temperature, ion mass 
fraction, and time after explosion are degenerate. 
We fix the excitation temperature to 10, 000K and 
the time after explosion to 1 day and vary only the 
ion mass fraction to represent the value of these 
combined parameters. The abundance scale fac¬ 
tors are limited such that the total optical depth 
stays below 10 4 to prevent unobserved Ca II fea¬ 
tures from appearing in the synthetic spectra, e.g. 
a A5020 doublet. 

We fit the model spectra to the obser ved spec¬ 
tra SN 2011fe from Parrent et al.l ( 2012i ) over the 
period —16.54 d to +2.87 d. For fitting purposes, 
normalization factors N s and N° are computed 
such that 


pX u p\ u 

/ N s T s x d\= / N°J T x d\ = 0.5, (1) 

A, 

where T\ is the specific flux density from either 
our synthetic spectra (Tp) or observed spectra 
(J 7 ^). We have chosen the bounds to be A/ = 
4500 A and A„ = 9000 A to cover the majority 
of the optical spectrum. A multi-variable simplex 
fitting routine is used to minimize the variance 
er 2 = ((N s jF x — N 0 jF x ) 2 ) over the range 7750A and 
8400 A, which covers most of the CaNIR absorp¬ 
tion feature. The CaNIR feature is chosen because 
it is the most conspicuous HVF, and it allows a fit 
without concern of blending with Si, as may oc¬ 


cur in the Ca H&K feature (IMaguire et al.l 12012 : 


Folevl 20131 : Silverman et alJl2015 1. Details of the 
dynamic models and line -fitti ng: pr ocedu re will be 
given in Mulligan fc Wheeler ( in preparation} ). 


4. Results and discussion 


Synthetic spectra for the models are shown in 
Figures [T]and [2] for the case of the supernova-only 
model (SN-O) and supernova with a shell of mass 
0.005 M 0 (SN+S), respectively. These figures also 
present the o b serve d spectra of SN 20life from 
Parrent et al.l 1 2012 + The photospheric velocity 
and temperature, abundance scaling factors ( C s 


and C E ), and variance of the fit (a) are listed in 
Table m The photosphere velocity and tempera¬ 
ture are based only on the fit to the CaNIR feature 
and are therefore not strong indicators of the ac¬ 
tual photospheric parameters; fitting of additional 
absorption features will give more reliable values. 
In the discussion b elow, phase is bas ed on Bmax 


on JD 2,455,814.4 (|Vinko et al.ll2012h . 


The SN+S model tends to produce a better fit 
to the width of the CaNIR feature over the en¬ 
tire period studied, although the variance usually 
differs by only about 0.1 between the SN-0 and 
SN+S models. The largest exception is for the 
earliest available spectrum (—16.54 d), when the 
SN+S model is a better fit by about 0.5 dex than 
the SN-0 model. This is also apparent by visu¬ 
ally inspecting the resulting spectra: the SN-0 
model leaves gaps both red-ward and blue-ward 
of the feature minimum. Part of the reason the 
earliest fit of the SN+S model is better than the 
remainder of the fits is the success in capturing the 
whole of the P Cygni peak; with the exception of 
-10.54d, there is excess flux at the peak in the ob¬ 
servations that the models do not reproduce. This 
is due to the choice of a fixed upper limit for fit¬ 
ting at 8400 A, which does capture the peak on 
—16.54 d, but is blue-ward of the peak at all later 
times. From visual inspection, the fit tends to re¬ 
sult in the flux near 8400 A matching the observed 
value but with an incorrect slope. A dynamically 
chosen upper wavelength limit for the fit would 
improve this result; this will be incorporated in 
future work. 


At the earliest times in the SN+S model, the 
photosphere lies in the shell rather than in the 
ejecta. The contact discontinuity between ejecta 
and shell material lies at 20,880 kms -1 . On 
—16.54 d our fit finds a photosphere velocity of 
21,660 kms -1 , which is within the shell; one day 
later the photosphere lies at 20,820 kms -1 , just 
inside the ejecta. This results in no measurable 
PVF in our model in these first two epochs. By 
— 13.21 d and later, the photosphere has receded 
into or below the freshly synthesized Ca, giving 
distinct PVF and HVF due to the ejecta and shell, 
respectively. At —7.27 d and later, the shell has 
become optically thin to CaNIR and the HVF due 
to the shell becomes very weak. 

There is a feature near 8000 A in both the ob¬ 
served and the SN+S model spectra that is miss- 
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ing from the SN-0 model spectra. This feature 
is due exclusively to absorption within the ejecta 
in the SN+S model. The lack of such a feature in 
the SN- O model suggests that it is a result of the 
interaction between the ejecta and shell. We have 
tested simplified models of a density enhancement 
at the reverse shock in the ejecta that do not im¬ 
mediately account for this interaction feature (IF). 
We will further explore the physical origin of the 
IF in our upco ming paper, Mulligan fc Wheelerl 
( in preparation! ). 


meas urable PVF, al t hough a PV F is identified by 


4.1. Velocity evolution 

Generating spectra with C s or C E set to a 
small value (10 -20 ) allows us to individually mea¬ 
sure the effect of the material in the shell and 
the ejecta. We flatten the ejecta-only and shell- 
only spectra using a continuum from the photo¬ 
sphere temperature and velocity, then identify the 
point of minimum flux of the CaNIR feature for 
the shell, which we refer to as the HVF veloc¬ 
ity, and ejecta, which we refer to as the PVF 
velocity. We report these velocities in Figure [3] 
with the SN 2 0 life P VF and HVF velociti e s from 
IParrent et al. ( 2012 ) and Silverman et al.1 ( 20151) 
shown for reference. 

The IF is identified by visual inspection of 
the ejecta-only spectra generated from the SN+S 
model. At times before —7.27 d, where there are 
multiple small features blue-ward of the CaNIR 
minimum, we chose the feature that had a wave¬ 
length range similar to the IF on the successive 
epoch. We find a local minimum in the flux at the 
IF and report that velocity in Figure [3] 

Comparing the HVF, P VF and IF ve l ocities 
in this work to those of Parrent et ah (2012) 


and ISilverman et al. ( 20151) requires recognition 
of the differing methods f or velocity measure¬ 
ment. Parrent et alJ (120121 ) report the v rnln pa¬ 
rameter they use for a SYNAPPS fit, which is 
equivalent to our photospheric velocity for the 
ejecta, but we have no equivalent for the shell 
oth er than the veloc i ty of t he contact discontinu¬ 
ity. Silverman et al. ( 2015 ) use a central velocity 
from a Gaussian profile line-fitting routine, which 
should be similar to our measured velocity, though 
some difference can be expected because the bot¬ 
tom of the CaNIR feature is not smooth. 

At the earliest times, the photosphere is in the 
shell so the SN+S model has only a HVF and no 


both 

-— - ' - 7 j " 

Parrent et al. 

(2015 

. The use 

Parrent et al. 

2012) 


(2012|) and ISilverman et al.l 
of the v m i n parameter by 


ence of a PVF that is absent due to very low 
effective optical depths. The method for iden¬ 


tifying PVF and HVF used by ISilverman et al 


(12015 ) and others may artificially require a two 
component fit because of the non-Gaussianity of 
the spectral feature as well as the choice of the 
P Cygni peak as the basis for the red-ward con¬ 
tinuum. The similarity of the velocities produced 
by the SYNAPPS and Gaussian-fitting methods at 
the early times cannot be discounted, however our 
model results call for deeper consideration of the 
techniques used to characterize the early data. 

The SN-0 model produces a PVF with a ve¬ 
locity t hat is higher than th e PV F velocity identi¬ 
fied b y IParrent et alJ ( 20121) and Silverman et al 
( 20151) by about 5,000 kms -1 at — 10.17d and ear¬ 


lier, but is too slow to explain the HVF. The PVF 
velocity of th e SN+ S mo del agrees with t hat o f 


Parrent et al.l ( 20121 ) and Silverman et al. ( 20151 ) 


on the first day that it is measurable (—13.21 d), 
but remains at about the same velocity over the 
following six days. At the la test times, it is again 
in v ery good agreement wit h Parrent et alJ ( 2012 ) 
an d ISilverman et, al.l (|2015h . 


T o explain the HVF in IParrent et al.l (|2012r ) 


and ISilverman et al. ( 2015 ), we suggest that the 
observationally-identified HVF is generated by 
the shell at the earliest times (—10.54d and earlier) 
and by the IF at later times. Our HVF velocity is 
about 3,000 km s -1 slower than their HVF veloc¬ 
ities during the shell-dominated phase; a slightly 
lower mass shell may explain the discrepancy, or it 
may be a result of the differing velocity measure¬ 
ment methods. We will further explore the effect 
of the mass of the shell in_our upcoming paper 
Mulligan fc Wheeler ( in preparation! ). 


4.2. Implications for progenitor models 

Our results disfavor progenitor models that pre¬ 
dict, or are consistent with, very little CSM. Mod¬ 
els in this category could be spin up / spin 
down models in which mass transfer has long since 
cease d (Di Stcfan o et al.ll2011tlDi Stefano fc Kilic 


201 2; Justham 20111 ) or models in which isolated 
white dwarfs explo de by pycnonuclear reactions 
( Ghiosi et al.! 201 5(1. In the absence of a CSM, 
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Table 1 


Parameters and variance of fit 


Model 

MJD a 

Photosphere 

velocity 

Phase b (lOOOkms -1 ) 

Photosphere 
temperature 
(1000 K) 

log C E 

log C s 

log a 


797.66 

-16.54 

24.80 

8.42 

1.98 


-1.73 


798.76 

-16.24 

22.26 

8.00 

1.14 


-1.60 


800.69 

-13.21 

17.97 

8.00 

0.43 


-1.89 


801.24 

-12.66 

19.58 

11.76 

0.06 


-2.14 

SN-0 

803.36 

-10.54 

19.85 

9.14 

-0.30 


-2.31 


803.73 

-10.17 

19.30 

8.00 

-0.81 


-1.93 


806.63 

-7.27 

11.95 

8.71 

-1.90 


-2.53 


811.73 

-2.17 

9.39 c 

19.55 c 

-2.17 


-3.16 


816.77 

2.87 

5.97 c 

22.52 c 

-2.12 


-2.76 


797.66 

-16.54 

21.66 

8.01 


0.41 

-2.32 


798.76 

-16.24 

20.82 

8.00 


0.08 

-1.65 


800.69 

-13.21 

18.35 

8.03 

0.12 

-0.27 

-1.72 


801.24 

-12.66 

16.98 

10.19 

-0.45 

-0.60 

-1.82 

SN-S 

803.36 

-10.54 

18.59 

8.18 

0.09 

-0.91 

-2.21 


803.73 

-10.17 

17.88 

8.00 

-1.33 

-1.63 

-1.81 


806.63 

-7.27 

19.09 

8.47 

-0.29 

-2.72 

-2.43 


811.73 

-2.17 

7.34 c 

18.90 c 

-2.39 

-3.05 

-3.19 


816.77 

2.87 

7.97 c 

22.81 c 

-2.18 

-4.11 

-2.79 


a Relative to JD 2,455,000.5. 

b Relative to B-band maximum on MJD 813.90 ( Vinko et al.ll2012h . 


c The photosphere is at a velocity that is interior to all of the freshly-synthesized Ca 
and thus the velocity is not well constrained by this fit. The high temperature is a 
result of the low velocity. 
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these models would not produce HVF by the 
mechanism modeled here. 


An important category of SN la explosion mod¬ 
els involve edge-lit double detonations whereby 
the explosion is triggered in a sub-Chandrasekhar 
white dwarf by the detonation of a thin layer of 
He on the outside of a C/O core. The required 
mass of He is estimated to range from ~ 0.05 to ~ 
0.1 Mq ( Fink et all 201 (1 fWooslev fe Kasenll201 ill , 
although the mass of He might be consi derab ly less 
if the He is enriched with C/O ( Shen fe Moorel 
12014 ). In the context of the current models, we 
note that the high-velocity ejection of such a He 
shell is unlikely to produce an observed HVF be¬ 
cause the velocity will be too high and the optical 
depth too low. 

Some progenitor model components can also be 
discounted as sources of the CSS. A thin accre¬ 
tion disk with a small solid angle would not ac¬ 
count for the ubiquity of the HVF and such disks 
are typically much less massive than the shell that 
we have invoked. A steady-state wind would re- 
sult in excess optica l light, which is not observed 


( Gerardv et al. 20041) . and would contain substan¬ 


tially less mass than does the shell presented here. 
The CSS may be a result of a thick disk or other 
phenomenon. Investigation into the composition 
of the CSS will provide insight into these sources. 


5. Conclusion 


We have compared the synthetic spectra from 
a model with only supernova ejecta (SN—O) and a 
model with the supernova ejecta colliding with a 
circumstellar shell of mass 0.005 -Mq (SN+S) to the 
observed Ca II ne ar-infrared triplet ( CaNIR) fea¬ 
ture in SN 2011fe ( Parrent et al.ll2012l) . The SN+S 
model generally produces a better fit to the width 
of the feature than the SN-0 model. The SN+S 
model spectra and observations of SN 2011fe have 
a sub-feature near 8000 A at late times which is 
absent in the SN-0 model spectra. This feature, 
which we dub the interaction feature (IF), is a re¬ 
sult of the interaction between the ejecta and the 
shell in our model. 

The combination of the shell-HVF and IF in 
the SN+S model more consistently reproduces the 
observed HVF velocity evolution than does the 
SN-0 model. The PVF velo city evolution in the 
SN+S model matches that of Parrent et al. ( 2012 ) 


and Silverman et al. ( 20151) when it first appears 
in this model and again near Brnax, but is about 
2,000 km s _1 fast at intermediate times. The 
SN-0 model fails to account for either the HVF 
or the PVF velocity at any time before about —8d, 
and accounts for only the PVF velocity after this 
epoch. 


We have made a plausible case that a shell of 
mass ~ 0.005 Mq accounts for the common pres¬ 
ence of HVF in SN la. Discrepancies between the 
spectra generated by our models and the observed 
SN 2011fe features may be improved by using a 
more finely tuned mass or structure of the shell. 
Fitting of additional absorption features is neces¬ 
sary to place constraints on the physical quantities 
of the shell, such as the geometry, mass, and com¬ 
position. These quantities will give us insight into 
possible origins of the shell and the cause of Type 
la supernovae. 


Special thanks to Jeffrey Silverman, Howie 
Marion, and Jozsef Vinko for many discussions 
about high-velocity features and the fits shown in 
this paper. Thanks to Vadim Gamezo for provid¬ 
ing the source model of the post-explosion super¬ 
nova. This work was supported in part by NSF 
grant AST-1109801. The Texas Advanced Com¬ 
puting Center (TACC) at the University of Texas 
at Austin has provided HPC resources that have 
contributed to the research results reported within 
this letter. 
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Wavelength [A] 
(a) -16.54 d 


Wavelength [A] 

(b) -15.24 d 


Wavelength [A] 

(c) -13.21 d 





(d) -12.66 d 


(e) -10.54 d 


(f) -10.17 d 



(g) -7.27 d 



(h) -2.17 d 



(i) +2.87 d 


Fig. 1.— Results_of fitting t he Ca NIR feature in the SN-0 model to SN 2011fe. SN 2011fe spectra are 
shown in black (Parrent et ah 2012 ) and the generated spectrum from the model is shown in red. The gray 
shaded regions at —13.21 d and later represent the approximate range of the IF in the SN+S model spectra; 
there is no equivalent feature in the SN-0 model spectra in these regions. 





















(a) -16.54 d 




(c) -13.21 d 





(d) -12.66 d 


(e) -10.54 d 


(f) -10.17 d 





(g) -7.27 d 


(h) -2.17 d 


(i) +2.87 d 


Fig. 2.— Resul ts of fitt ing the CaNIR feature in the SN+S model to SN 2011fe. SN 2011fe spectra are shown 
in black ( Parrent et alJl2012h . the generated spectrum from the model is shown in solid red, a spectrum 
showing the effects of only the shell is shown in dashed blue, and that of only the supernova ejecta is shown 
in dashed red. At early and late times the ejecta and shell component of the spectrum, respectively, is shown 
as a continuum due to high photospheric velocity (early times) and thinning of the material (late times). 
The gray shaded regions at —13.21 d and later represent the approximate range of the IF in the SN+S model 
spectra. 
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Fig. 3.— Evolution of the Ca II near IR triplet feature velocities for the SN-0 model (blue, dash-dot, 
■), the SN+S model (das hed, PVF: red , ▼; H VF: yellow, ▲; IF: magenta, ♦), a nd the Ca II PVF and 
HVF velocities reported bv IParrent, et, al. ( 20121) (x. PV F: green; HVF: brown) and Silverman i et_aL ( 2015 1 


(+, PVF: cyan; HVF: dark brown). The Parrent, et al.l data stop at —10 d, but the Silverman et al. HVF 


velocities extend to Bmax where they correspond to the IF velocity in the SN+S model. 
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